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Abstract 

We study the dynamics of entanglement entropy for weakly excited states in conformal field 
theories by using the AdS/CFT. This is aimed at a first step to find a counterpart of Einstein 
equation in the CFT language. In particular, we point out that the entanglement entropy 
satisfies differential equations which directly correspond to the Einstein equation in several 
setups of AdS/CFT. Wc also define a quantity called entanglement density in higher dimensional 
field theories and study its dynamical property for weakly excited states in conformal field 
theories. 



1 Introduction 

The AdS/CFT correspondence [U [2] is a remarkable fundamental relation which connects gravita- 
tional systems and quantum field theories as an equivalence. In spite of its recent successful de- 
velopments, we are still far from the complete understanding of the basic mechanism of AdS/CFT 
correspondence (or gauge/gravity duality). The aim of the present paper is to report a modest 
progress in this direction. In particular we would like to study what the Einstein equation in the 
gravity side corresponds to in the quantum field theory side. 

Since the AdS/CFT correspondence relates gauge invariant quantities between both sides, the 
Einstein equation itself, which is written in term of the spacetime metric, is not directly interpreted 
in the dual quantum field theory. Therefore we need to find a counterpart of Einstein equation for 
gauge invariant quantities. We argue that the holographic entanglement entropy (HEE) is one of 
the best quantities for this purpose. 

The entanglement entropy in quantum field theories and more generally quantum many-body 
systems has been intensively studied recently (see e.g. the review articles O IU O [6] ) . In AdS/CFT 
we can holographically calculate the entanglement entropy in a gravity dual as an area of min- 
imal surface as conjectured in [71 [8]. This holographic entanglement entropy (HEE) calculation 
was proved in [9] quite recently by using the bulk to boundary relation [2]. See [10] for a proof 
when the subsystem is a round ball. Also refer to [111 [T2| [13] for strong supports towards a proof 
within AdS3/CFT2. This holographic calculation of entanglement entropy was also developed in 
time-dependent backgrounds [H] , which has been applied to the quantum quenches [El [HI [17] , 
the de-Sitter space [20] and energy flow [21]. In [16], a falling particle in AdS was considered as 
the holographic dual of local quenches [19] and its HEE has been computed. In [T7], an analyt- 
ical framework for holographic counterpart of global quantum quenches and their entanglement 
entropies in CFTs [18] has been discovered. In this paper we will study how a small perturbation 
of HEE evolves dynamically by solving the Einstein equation in AdS spaces. A final goal will be to 
rewrite the perturbative Einstein equation in terms of the HEE and we will do this explicitly in sev- 
eral examples of AdS3/CFT2. At the same time, our results describe the behavior of entanglement 
entropy for weakly excited states. 

Note that the complete information of HEE for arbitrary subsystems is essentially the same as 
that of the spacetime metric in the gravity dual (see e.g. [23]). For this correspondence we do not 
need to know the details of the Lagrangian of matter fields coupled to the Einstein gravity. On the 
other hand, if we want to reproduce the spacetime metric only from the information of holographic 
energy stress tensor [22] we need to employ the precise form of Einstein equation and thus this 
requires the details of matter fields. 



When the size of subsystem is small, we can find a simple relation between its entanglement 
entropy and the total energy inside it for excited states. This is called the first law-like relation 
and has been first obtained in [25] when excited states are static and translationally invariant. Our 
analysis in this paper provides a proof of this first law relation for spherical subsystems in the 
presence of time-dependent excitations. 

We will also study the dynamics of entanglement density introduced in [16] for two dimensional 
field theories. We will extend this quantity in higher dimensions. Then we consider evolutions of 
entanglement density imposed by the Einstein equation. 

This paper is organized as follows: In section two we will present a general strategy of conducting 
our perturbative analysis of HEE and review the first law-like relation. In section three, four and 
five, we present our analysis of HEE in AdS3/CFT2, AdS4/CFT3 and AdS5/CFT4, respectively. In 
section six, we study the higher dimensional entanglement density and its perturbations. In section 
seven, we summarize our conclusions. 

2 Perturbative Calculation of Holographic Entanglement Entropy 

2.1 General Strategy 

We consider a perturbation of the pure AdS^+i metric in the Fefferman-Graham (FG) gauge as 
follows: 

The coordinate x^ (^ = 0, 1, 2, • • •, d— 1) describes the d dimensional Lorentzian space i?^''^~^ where 
the d dimensional conformal field theory (CFT^^) lives in. The parameter K describes the radius of 
the AdS space. 

We consider the perturbative expansion of the form: 

9[iv = V/j-u + h^y, (2.2) 

assuming that h^y is very small. We are interested in only the linear order of h^j^y. The dynamics 
of this perturbation is determined by the Einstein equation as usual 



7^ab - ^nga, - ^^^^^^Qa, = T^,\ (2.3) 



(G) 
where Tfi^ is the energy stress tensor for the matter coupled to the Einstein gravity (a, 6 = 0, • • •, d 

are indices of coordinates of the d + 1 dimensional space). 

The holographic entanglement entropy (HEE) in a d -|- 1 dimensional AdS gravity is given 

in terms of the area of d — 1 dimensional extremal surface jj^ in a given spacetime so that the 



boundary of 7^ coincides with that of the subsystem A [TJ [H]. In this paper we are interested in 
the first order correction of HEE under a metric perturbation. This can be conveniently computed 
as follows. First we start with the extremal surface 7^ whose shape is already known in the pure 
AdS and calculate its area. Because the pure AdS spacetime is static, ja is a minimal area surface 
on a canonical time slice. Next we evaluate the area of the same surface "ja in the perturbed metric. 
The differen between these two gives precisely the first order correction of HEE AS a- Thus we do 
not need to know how the shape of the extremal surface is modified under the metric perturbation. 
This simply stems from the fact that ja satisfies the extremal surface condition in the pure AdS 
space. 

Therefore we can calculate the shift of HEE AS a due to the metric perturbation ()2.2p from the 
formula 

^Sa = ^ |(dC)'^^i v/gWg^gW"^ (2.4) 

where ^s a coordinate of the d — 1 dimensional extremal surface ja as employed in [16]. G^^' 
and G^^' represents the induced metric on ja with respect to the pure AdS and its first order 
perturbation. In this paper we consider examples where the subsystem A is given by a round ball 
for which we know the analytical expression of the surface ja in the pure AdS. 

2.2 First Law-like Relation 

A useful property which is enjoyed by the perturbed HEE is the first law-like relation [25]. It takes 
the following form: 

T,ff ■ ASa = AEa. (2.5) 

When we choose A to be a round ball with radius I, the entangling temperature Tf^ff takes the 
universal value 

T.„ = ^. (2.6) 

Also AEa denotes the total energy in the region A and is written as 



AEa = / Ttu (2.7) 

where Tu is the energy density in CFT^;. If we perform the expansion of h^i, defined in (j2.2|) as 

V = ^'^H^^ + • • •, (2.8) 

in the near AdS boundary limit z — )■ 0, the holography energy stress tensor can be obtained as 
follows [22] 



IQttGn 



V = -TT-TT-H,.. (2.9) 



This relation ()2.5|) was originally confirmed in static, isotropic and translationally invariant 
perturbations of the metric. Therefore it is intriguing to see if it holds without this assumptions. 
This is another motivation of this paper. 

3 Analysis of Perturbed HEE in AdSs/CFTa 

To study the HEE in AdS3/CFT2, we set {x'^,x^) = {t,x) in (j2.ip . We choose the subsystem A for 
the entanglement entropy Sa to be an interval —1/2 + S, < x < 1/2 + ^. Then the corresponding 
minimal surface 7^ is parameterized as 

. II , ^ 

t = const., X = 4 + -sin 99, z = -cosip. (3-1) 



Then the shift of HEE (j2.4p reads in terms of the perturbation of the metric h^y-. 

1 r^^ , fl/2 cos (^ . 2 , 2 , 

^^A = -Trpr- I dip —5 cos iphxx + sm iphzz -2smip cos tph, 

8Gn y_^/2 4:Z^ 

R r'^ 



dip cos (fhxx, (3.2) 

oLrN J-n/2 

where we employed the EG gauge in the final expression. 

3.1 HEE in AdSg Pure Gravity 

As the simplest example in AdS3/CET2, we would like to calculate AS a in the pure Einstein 
gravity for AdSs. The equations of motion for the metric perturbation h^u reads 

dz{dxhtx - dthxx) = 0, 

dzidxhtt - dthtx) = 0, 

{dz - zdl)hxx = {dz - zdl)htx = {dz - zdl)htt = 0, 

dz{htt - Kx) - zdlhtt + 2zdtdxhtx - zd'^h^x = 0. (3.3) 

Note that this system is topological in that there are no propagating degrees of freedom as usual 
in three dimensional pure gravity. 

By requiring that h^^{z,x) is order O(z^) so that only normalizable modes are excited, these 
equations can be solved as follows: 

htt = hxx = z^H{t,x), dfhtx = z'^dxH{t,x), dxhx = z'^dtH{t,x), (3.4) 

where the function H{t,x) satisfies 

id^-dl)H{t,x)=0. (3.5) 



Finally, by using this solution, (|3.2p is rewritten as follows: 



T>l2 j-TT/Z / I \ 

ASAi^,l,t) = -—r- / dvcos^<fH[t,^ + -sm^) 



-■it/2 
72 /"7r/2 / 7 

— - / difcos^ ip-Tu[t,i + -smtp] , (3.6) 



4 7_V2 - "^- 2 

where we employed the relation H{t,x) = ^j^^ Ttt{t,x) by setting d = 2 in (12. 9p . 
By taking the Fourier transformation 

/oo 
dxe-'''''H{t,x), (3.7) 

-oo 

we obtain 

i? 2sinf -Z/ccosf , , , , 

A5a(A;, i. = ^^ ^^3^ ■ m, t)- (3.8) 

In this way, we find that ASa is related to the metric via this non-local transformation. Moreover, 
it is straightforward to see that it satisfies 



«? + (t-| 



ASAik,l,t) = 0. (3.9) 



In summary, we find from (|3.5|) and (j3.9p that AS a satisfies the following "equations of motion 
for entanglement entropy": 

id^-dl)ASAiCJ,t) = 0, (3.10) 



ASA{^J,t) = 0. (3.11) 






The first equation (j3.10p shows that the quantum entanglement propagates at the speed of light 
in the x direction. The second one (j3.1ip describes an evolution in the width (or radial) direction, 
which is analogous to the wave equation in an AdS spacetime. We believe that the presence of two 
constraint equations for AS a as in (j3.10p and (j3.1ip is peculiar to the AdS3/CFT2 duality and this 
is due to the fact that the gravity does not have propagating degrees in three dimension. 

3.2 HEE from Einstein-Scalar Theory 

Since the pure AdSs does not have any propagating degrees of freedom, it is more interesting to 
consider a Einstein-matter theory on AdSs. Though we will not write down all components of 
Einstein equation, we would like to note that the tt component reads 

d,h,. - zdlh,, - 2zTP = 0, (3.12) 

(G) 
where T^^, is the energy momentum tensor in the d + 1 dimensional gravity. This should be dis- 
tinguished from the holographic energy tensor T^u for the d dimensional dual CFT. By integrating 
this equation, the shift of holographic entanglement entropy ASa (13. 2p is expressed as follows: 



A5A(i, C, = -^7:^ / dipR cos iph^^ f - cos y?, t, ^ + - sin if j 

= -^/ dy^cos^ -2 dz" dz'r,f z',t,e + ^sin(^ — (3.13) 

oCzTv j_| L io Jo V Z J Z' 

P / I \ 2 ' 

+ ^H \t,^ + -sm(p j COS If 

In particular, we consider the Einstein-scalar theory which is defined by a free scalar field (p 
(mass m), which is minimally coupled to the Einstein gravity. The Einstein equation is given by 
()2.3p with the energy stress tensor is given by 

^/i^^ = ld^<t>du4> - ^g^u [{d^? + m^cp'] , (3.14) 

where we normalize the scalar field appropriately. 

We consider a perturbation of the scalar and AdSa metric in the EG gauge as follows: 

(j){z,t,x) = yfei){z,t,x), 

(3.15) 
g^u{z,t,x) =r]^^ + €hfj,^{z,t,x), 

where e is an infinitesimally small parameter of the perturbation. 
The equation of motion of this perturbation is given by 

z (z (V(°'°'') {z,t,x)- V(°'''°) {z,t,x)+ V'^''^'"^ {z,t,x))- ^(^'O'O) {z,t,x))- m'R^i;iz, t,x)=0. 

(3.16) 
After the Fourier transformation, the normalizable solution for this equation is given by 

/oo poo p—itui+ixk , , 

du: / dk-- 1. , a,l^ (<^(^' ^)) Jv ra^RUl (^V^^ " kA , (3.17) 

-oo J —oo (oj^ — k )^ 

where {0{uj, k)) denotes the expectation value of an operator dual to the scalar field (j) (up to a 
normalization factor) ptl24|. 

The entanglement entropy is decomposed into two parts: 

ASAiL I, t) = A5^(C, /, t) + A5f (C, /, t). (3.18) 

The AS^ is the same as that in the pure AdSa case and is given by (|3.6p . The other term AS"^ 
comes from the contribution from the matter field and is expressed as 

77> r~ f cos <~p rz^' / i \ // 

AS^{^,l,t) = -^j[ d^JJ dz" j^ dz'rP (^z',e + ^sin^,tj . '- 
— I dip cos If / dz" / dz' / duiidoj2dk\dk2 



I6G7V 

X zV (0(wi, fci)) (0(^2, A;2)) F(^',c^i,a;2, A;i, A:2)e<"('^i+'^^)*+('=i+^2)(?+l--^)) 



F{z,uji,u}2,ki,k2) is given by 

F{z,UJi,0J2,h,k2) = — {ul - kl)^ [ul - kl)^ 



'^\l'^\ - k\jy^i ( 2:^11^1 - ^1 ) ( (1 - J^) «/;/ f -Zy ^2 ~ ^2 ) + ^\l'^ 



-J^ 



LO2 — /cl'/jy-l ( Z\l UJ2 — k^ 

'y I z\lij\ — k\\ \Jv\ zy ^2 ~ ^2 ) ■ {—kik2z'^ + 2^"^ — 2v — ujiuj2z'^) 



Z\ IJJ2 — ^2 (^ ~ 1) ' Jy-l I -2\/ ^2 — k 



(3.20) 
(3.21) 



where we defined v = \/m?B? + 1. 

We can also act a differential operator to get rid of the AS'^ contribution and obtain the 
following constraint equations which are satisfied by AS*^: 



(a2-a|)A5A(/,t,c) 

duJidU2dkidk2F^^\ki,k2,UJi,L02, I) {0{uJi,ki)) {O{l02, ^2)) e^'^^^+^^^'+'^''^+''^^^ . 



(3.22) 



df-\df-^]ASA{l,t,0 



I 



da;idL^2dA;idfc2F(2)(fci,A;2,^i,W2,0 (OK,^i)) (0(W2,A;2)> e-^("i+"2)*+'('=i+'^2)«, 



(3.23) 



where we did not write explicitly F^^' and F^"^' as they can be easily obtained from (j3.19p . 

In summary, we have obtained a counterpart of the perturbative Einstein equation in terms 
of entanglement entropy. Therefore we can regard ()3.22p and (|3.23p as the perturbative equations 
of motion for the entanglement entropy. Once we specify the expectation value of energy density 
Ttt{t, x) and the scalar operator 0{t, x) as functions of t and x, then the differential equations p.22p 
and (j3.23p determine the time evolution of A5^(i^,/,t). 

They share a similar structure with that of the Einstein equation because the left hand side 
comes from the gravitational physics (-H- entanglement) and the right hand side corresponds to the 
matter fields (-H- operator expectation values). Note also that they are manifestly gauge invariant 
(as far as we fix the AdS boundary coordinate), as opposed to the Einstein equation where the 
metric changes under the coordinate transformation. 

One may notice that we only used the Einstein equation involving the space-like component 
hxx- In order to take into account the one for the time-like component we need to consider the 
entanglement entropy for boosted subsystems. We will leave the details of these for a future 
problem. 



3.3 Proof of First Law-like Relation 

When I is very small, we can neglect the contributions from the matter fields ^S^ . Therefore in 
this limit, we find 

TJj2 

ASA{tl,t)c^^^Hit,C). (3.24) 

We can show this relation (|3.24|) confirms the time-dependent version of the first law-like relation 
(12. 5p of the entanglement entropy for AdS3/CFT2. Indeed, since by definition we find tu = H{t, x) 
we can easily reproduce (|2.5|) from (|3.24|r . 

4 Analysis of Perturbed HEE in AdS4/CFT3 

Now we move on to higher dimensional cases and here we especially consider a AdS4/CFT3 example. 
Since the pure Einstein gravity is already dynamical in four or higher dimension, we will concentrate 
on the pure gravity just for simplicity. 

4.1 Einstein Equation 

We consider a metric perturbation of the AdS4 space. The metric is again given by (j2.1|) and (j2.2p 
with {x^,x^,x'^) = (t,x,y). We require that h^^, is order O(z^) in the limit z — )• so that we can 
keep only normalizable deformations. 

By performing the Fourier transformation 

hf,,{t,x,y,z) = ^ I dudk^dkye-"^'+'''=^''+"'yy ■ h^,{z,LO,k^,ky), (4.1) 

the perturbative Einstein equation leads to 

za2 V - 25, V + (^^ - ^x - ^^)^V = 0> (4-2) 

for all components of perturbations hu, htx, hfy, h^x, hxy, hyy. This is easily solved as 



h^u{z,u,kx,ky) = 'iJ- ■ H^y{u:,kx,ky)- {u? -k^) * • 2:2 • J^i^[y/ ufl -k'^z). (4.3) 



Note that the explicit form of the Bessel function reads 




^3/2(x) = V;;^f-cosx + ^). (4.4) 



vrrr y x 

We can show the following behaviour near the AdS boundary z — )• 0: 

.3 



hf^:y{z,u},kx,ky) c^ z ■ Hf,„{u},kx,ky). (4.5) 



^note that due to our conventions for the AdSz/CFT2 case we should replace in (|2.6p / — i- 1/2. 



Since we are interested in non-singular and normalizable solutions, we restrict to the range w^ > fc^. 
We can show that the Einstein equation is equivalent to (14. 2 p and ()4.6p 

** ~ T2 7^2 [~'^^xkyhxy + (Kj. — ky)hyy\ , 
htx = ^,2 _ fc2^ [-hi^"^ + kl)hxy + kx{u}^ - kl)hyy] , 
illy (^ f^x^xy y yy}"! 
^^^ = I2 2 [-'^kxkyhxy + {u? - ky)hyy] . (4.6) 



In particular, we can find the relation from ()4.6p 

ht = hxx + hyy, (4.7) 

which is interpreted as the traceless condition of energy stress tensor in the dual CFT. 

4.2 Calculations of HEE 

We consider the shift of HEE AS a by choosing the subsystem ^ to be a disk with a radius I. In 
the pure AdS4, the HEE is computed as the area of the minimal surface ja given by the half of 
sphere parameterized by 

X = lsm9cos(j) + X, y = I sin 6 sin (j) -\- Y, z = lcos6, (4-8) 

with the range < ^ < 7r/2 and < (j) <27t. 

The shifted amount of the HEE, denoted as ASa due to the linear perturbation of the metric 
is found by using (j2.4p : 



ASAit,X,Y,l) 

n2 r2TT /'i"/2 ■ a 

= -— — / dcj) dO 2^ [(1 — sin 9 cos 4>)hxx — 2hxysin 6 cos (f) sin cf) + (1 — sin 6* sin 4')hyy'j . 

oGjsf Jq Jq cos 6 

(4.9) 

We take the Fourier transformation of ASA{t,X,Y,l) with respect to t,X,Y, which is denoted by 

ASAi^,kx,ky,l): 

ASA(.t,X,Y,l) = -i^ f dudkxdkye-'^'+"'-''+"'y'' ,ASA{oj,kx,ky,l) 
ASAiuj,kx,ky,l)= I dtdXdYe''^*~'^-^"'''y^ ASA{t,X,Y,l). (4.10) 



Then we find 



3R^ RF r^'^ ,, sin ( 
do K 



2-K 



8Gn V 2 Jo cos- u jQ 

I cos 9 \^/2 



dcbliS 9) • e«'s™e(fciCos0+A:ySin<7i) 



Voj'^ — k'^ 



■ J3/2{lcos9\/oo'^ — k'^), 



(4.11) 



where we defined 



(1 — sin 0COS (t))Hxx{^,k) — 2sm 9 cos (j) sin (pHxyiu), k) + {1 — sin Osin (j))Hyy{u!,k). 

(4.12) 

(4.13) 



We can perform (j) integral by using the formula of Bessel functions 



Jniz) 



1 

2^ 



27r 



^mm—izsmi 



as follows: 

1 

2^ 



2tt 



,0)-e 



il sin 9{kx cos (p+ky sin <; 



{HxX + Hyy) I 1 



sm 



ii^. 



Jo{klsin9) ^ sin^ (e^*" J„2(A:/sine) - e-^'"" J2{kl sin 9)) 



sm 



-{Hyy - Hxx) {e^'"J^2{klsin9) + e-'^''^J2{klsin6)) 



1 



sm 



Jo{klsin9) + ^^^ (a;2 - — ) J2{klsine) 



Htt. 



(4.14) 



Here the angle a was introduced such that kx cos(j) + ky sin<j) = —ksin((p + a) i.e. sin a = —kx/k 
and cosQ = —ky/k; we employed ()4.6p to get the final equation. 

Thus we can express ASa{^, kx, ky, I) in term of the (holographic) energy stress tensor Tu{uj, kx, ky) 
as follows 



ASA{uj,kx,ky,l) 

2x/2vr5/^/3/2 W2 



sml 



where Q{1) is defined by 



Q(0) = J3/2(/cos0\/tj2-A:2) 



sm 



Jo{klsin9) + 



sm 



/fc2 



VCOS^ 



a; 



Q(^), 



J2{klsine) 



(4.15) 



. (4.16) 



4.3 First law-like Relation and Translationally Invariant Limit 

By taking the limit / — t- in (j4.9p . we obtain 



„oZW /-vr/z 

ASa = — :^— / d6lsin6'cos6l(l-sin2 6'/2)(iJ^^ + i:fj/y), 



37ri?2/3 



32GAr 



i^*. 



(4.17) 
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where we employed the relation (|4.7p . On the other hand the total energy in A reads 

AEa = TTl^Tu = ^^Hu, (4.18) 

by using (j2.9p . Therefore we can confirm the first law-like relation (j2.5p with Tgjj = ;^. In other 
words the energy density and ASa are related to each other via 

ASa = -^Tu. (4.19) 

Moreover, it is also intriguing to note that we can obtain the same result (|4.19p when we take 
the translationally invariant limit kx, fcy — t- even if we keep / and oj finite. This fact can be shown 
by explicitly evaluating ASa as follows: 

ASA{uj,x,y,l) 

SttR^I^^ [^ r/\^ sine / .2^ /2^2gi^4^X 

4.A Metric Shift in CFTs 

In the main part of this paper we have considered only normalizable perturbations which correspond 
to the dynamical change of state in a given CFT. However, we would like to study non-normalizable 
modes only in this subsection. Indeed, when we consider the time-dependent background corre- 
sponding to quantum quenches induced by a sudden change of the metric, we need to take into 
account the non-normalizable modeo 

For this, as a solution to the Einstein equation, we assume 

h,,{t,x,y,z) = i^l . I" ^e-*^*^^,(a;)a;3/2^3/2^{i) (^^)^ (421) 

where the Hankel function is defined by 

TTZ \ Z 



^3/2(^) = \/-(--l)e*'- (4-22) 



Near the AdS boundary z = 0, it approaches 

'■^ doj 



h^,{t, X, y, z = 0) = / -^e-'^'H^.iuj). (4.23) 



^If we consider quantum quenches induced by excitations of matter fields such as scalar fields, we will have milder 
metric backreactions such as the Vaidya metric. Refer to e.g. [26) for recent developments of holographic quantum 
quenches. 
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Then the HEE behaves hke 

as'aCw) = p / de—= 2 - sm^ e ^ — myjiuj cos 6i)(f^^ + h. 



8V2GJ, Jo ""^^^^V 4 y ^^3/2V^--^n^^.x - ^^..^ 

= |§^ (^^ + ^^ V) (^..(a;) + Hyyiu:)), (4.24) 

where 6 is the UV cut off and is related to the lattice spacing a via 6 = a/ 1. 
The divergent part of the HEE is simply obtained as 

ASl»(t) = ^(4 - l'd^){H,.{t) + Hyyit)) ■ !-, (4.25) 

where a is the UV cut off (lattice spacing). This agrees with the expectation from the area law. 
The more non-trivial contribution is the finite term, evaluated as follows: 



-^finite/ 



A5:^™^^(i) = —— ■ df H^t, X, y) + Hyyit, x,y)]. (4.26) 



^^^ ■ Clt \^Ha,x[t, X, y) + t.yyy 

5 Analysis of Perturbed HEE in AdS5/CFT4 

Here we analyze the perturbed HEE in the AdS5/CFT4 setup. Since the calculations are parallel 
with the previous section, our presentation will be brief. 

5.1 Solutions to Perturbative Einstein Equation 

The Einstein equations for the Fourier transformed metric perturbations 

are equivalent to 

zdlh^y - ^dzh^y + [u? - kl - kl - kl)zh^i, = 0. (5.1) 

together with 



^« = ( 1.2^ {2kik2hi2 + 2kik3hi3 + {-kl + k1)h22 + 2k2k3h23 + {kl-kl)h33), 

_ -{klk2 + k20j'^)hi2 - {klk-i + ksu}^)hi3 + {kiu'^ - kik^)h22 - 2kik2k3h23 + {kioj"^ - /ci/c|)/i33 
*' ~ u;{-kf+u^) ' 

ht2 = {klhi2 + k2h22 + k3h23), 

hfs = ikihi3 + k2h23 + k3h33), 

CO 
^11 = / 1.2^ {2kik2hi2 + 2kik3hi3 + {k^-U^)h22 + 2k2k3h23 + {kl-UJ^)h33). (5.2) 
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The normalizable solution of the first equation (|5.ip is 

h^y{z,uj, k) = S{u? - kY^^^MzVoJ^ - k^)H^^{uo, k) = z^H^^iuj, k) + 0{z^) (5.3) 

5.2 Calculation of perturbed HEE 

Consider the HEE Sa for the subsystem A whose boundary is defined by a two dimensional round 
sphere S"^ with radius I. The corresponding minimal surface in the pure AdSs is parameterized by 

z{9, (j), p) = I cos 6, xi{9, (j), p) = IsmOcoscj) + Xi, 

X2(6',(^, p) = / sin sin cos p + X2, X3(0, <^, p) = / sin sin </> sin p + X3. (5-4) 

Then the variation of the entanglement entropy is computed as 



A5a 

i?3 



^/2 /■2^ sin^^sinc 



dcj) dO dp 
SGat Jq Jq Jq cos'^ 



ill 



sm 



cos (f>\ — 2sm(j)cos{(j)){hi2COsp + hi^sin^p)) 



-2/i23 cospsinpsin + /i22 



1 



sm 



cos p sin (/> 1 + /i33 



1 



sm 



2 - sin (/)sin p 



(5.5) 



After the Fourier transformation 

AS{uj,k,l) = f dtd^X AS{t,X,l)e''"^-''''^' 
we obtain 

AS{uj,k,l) 



m^ 



■K j-n/2 1-2-K 



^ / 9 lox I "'•r I dS I dp n — J2{lcos9\/bj'^ - A;2)x 

Gn{w^ -k^) Jq Jo Jo cos 6 

il sin 0{ki cos (f>+k2 sin (f> cos p+k'i sin (/> sin p) p^ U A q] 



(5.6) 



with 



P{u,kA,0,p) = Hu{^,k){- 



sm 



cos^ (j)) — 2 sin (f> cos (j){Hi2{uj , k) cos p + Hi^{ijj, k) sin(p))- 



- 27^23(0;, A;) cos p sin p sin (j) + H22{i^Tk){^—^ cos psin (j)) + H^^{oj,k){^—2 sin (psvn p) 



sm 



sm 



5.3 First Law-like Relation 

Using ()5.5p in the limit of small / (j2.8p . we can perform the integrals as: 



ASa 



15G7V 



d(j) de sin^ 6* sin cos 6l[-2(Fii + H22 + ^^33) + 2Fii sin^ 9 cos^ <^ 
io Jo 

■3) sin^ (^sin^ 0] 
{Hu+H22 + Hs-i). (5.7) 



+ (i:r22 + i:^33)sin2,^sin2 6l] 
2/4i?37r 
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Since from equation (|5.2p it is immediate to obtain Hf^ = Hn + H22 + -^33, we can finally show 
the relation : 

9/4 n3„ 

This proves the first law- like equation (j2.5p by using (j2.9p . 

More generally, we can show the relation (j5.8p if we take ^1,2,3 — ^ limit with keeping uj and I 
finite as in the AdS4 case. 

6 Entanglement Density and Its Perturbation 

Since the entanglement entropy measures the entanglement between a certain region and its out- 
side, it is a highly non-local quantity. When we are interested in local physics, the entanglement 
entropy sometimes smears its essential effect. To improve this, the quantity called entanglement 
density was introduced in [16] for two dimensional field theories. This quantity measures the en- 
tanglement between infinitesimally small two regions and is guaranteed to be positive due to the 
strong subadditivity relation. We will generalize this quantity in higher dimensional field theories 
and study the behavior of entanglement density in excited states by using the perturbative Einstein 
equation in the holographic description. 

6.1 Entanglement Density in Two Dimension 

Let us start with the entanglement density in two dimensions introduced in [16]. This is defined by 

n{l,i,t) = \{^^dl-df^SA{l,i,t), (6.1) 

where we chose the subsystem A to be an interval such that its width is I and its center is at x = ^ 
as in the holographic description (j3.ip . We also write the shift of entanglement density compared 
to that of the ground state as An. 

In our holographic setup of the pure gravity on AdS3, this is evaluated as follows 

A /, ,N R -2sinf + /A;cosf ^^^, ^ 

Thus in this case we find the following simple relation is satisfied for the small excitations around 
the ground state: 

A5A(e,t,/) + ^^-An(^,t,0 = 0. (6.3) 
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In the gravity dual of the Einstein scalar theory, we can show from (|3.23|) that this relation is 
modified: 

The first law-like relation (j3.24p can be expressed in terms of the entanglement density as 

limAn(C,t,/) = -^r«(C,t), (6.4) 

which can also be seen easily by taking I — t- limit in (16. 2p . In this way, the entanglement density 
is equivalent to the energy density when its width / is vanishing. Though this was already noted in 
|16] when Ttt{S,, t) does not depend on ^ and t, here we gave a proof of this for more general cases. 



6.2 Entanglement Density in Higher Dimensions 

Now we would like to define the entanglement density for field theories in higher dimensions. 
Consider a d(> 3) dimensional Minkowski spacetime and define its coordinate {t,x) (if is a d — 1 
dimensional vector). We specify an arbitrarily given subsystem ^ by x = x{C), where C" (o = 
1, 2, • • •, d — 2) is the coordinate of the boundary dA of the subsystem A. Remember that A is 
defined on a certain time slice t =const. The unit normal vector at a point on dA (toward the 
outside direction) is denoted by A^(C)- 
We deform the subsystem A as 

x = x{C) + 6xn{C)-N{C). (6.5) 

In this setup we define the entanglement density by 

''^^^^^'^ = -sxJc)'Ucr ^'-'^ 

assuming C ^ (^'. Notice that this quantity depends not only on (^ and C but also on the choice 
of subsystem A itself. In this sense, the entanglement density in higher dimensions {d > 2) is not 
completely a local quantity. 

We can show that this quantity is non-negative due to the strong subadditivity of the entan- 
glement entropy [271 128]. To see this, let us assume 6xn = Sxn + Sxn , such that only one of fe„ 
and 6xn can be non-zero (or equally 6xn • Sxn = 0) for any points on dA. Then the strong 
subadditivity leads to 

S{x + 5x^^'^) + 5(f + (5x(,2)) > 5(f) + S(f + feW + 54^)), (6.7) 

where S{x + 6xn ) denotes Sa for the subsystem A defined by x = x{() + 6xn {() ■ N{C)- By 
performing the Taylor expansion of ()6.7p . we find that n^(C, CO is positive. 
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6.3 Holographic Calculation in Pure AdS 

Now we present a holographic calculation of entanglement density in any dimension. We take the 
subsystem ^ to be a half space defined by y{= x^-i) < 0. We perform an infinitesimal deformation 
of A (|6.5p by choosing (^" = x" (a = 1, 2, • • •, d — 2). In this case, the normal vector N{Q is in the 
y direction and thus the boundary dA of the subsystem A is simply expressed as 

y = fe„(x"). (6.8) 

We would like to calculate Sa for this deformed subsystem by using the AdS/CFT. Consider 
the pure AdS^+i in the Poincare coordinate {z,t,xi, ,X(i-2,y)^ whose metric is written as 

ds^ = n^dz^-dt' + ^tld^l + dy\ (6.9) 

z^ 



Before we do the deformation ()6.8p . the minimal surface responsible for 5^ is simply given by 
the plane y = in (16. 9p . Therefore, we can describe the minimal surface after the deformation 
dEHD as 

2/ = /(z,f), (6.10) 

where x = (xi, a;2, • • •, Xd-2)- 

To calculate the HEE , we take a variation of the area functional in order to find the minimal 
surface: 



S 



4Gjv J z'^-i \ 
Up to the second order expansion of /, Sa is approximated as 



-(S)'^g(i:)^- 






By taking a variation with respect to /, we find the following equation 



fAM)-'^-''^ <-^' 



By performing the Fourier transformation of f{z,x), we obtain 

We can find solutions to this equation that do not diverge as we take z to infinity: 

f{z) = ziKa{kz). (6.15) 
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Also we can find the solution that satisfies /(O, x) = ^(x) at the AdS boundary is given by 
,, ^, T(d-l) f , , o C{x)z'^ 



7r^T{l)J ((x-x)2 + z2)rf-i 

Now let us substitute this to the area functional (j6.12p . We ignore the divergent part because 
its dependence of .^(x) is local and we are not interested in it. Therefore we concentrate on the 
finite part given by 

This is evaluated as follows: 

^ 87r^r(f )G^ J ' ' (^1 - ^2)2('^-i) 

In this way, we obtain the entanglement density from the AdS/CFT: 

6'Sa 



n{xi,x'2)j 



6ax-'i)6C{x-'2) 
R^-^d-Vi^d-l) 1 



(6.19) 



4vr%^r(f )G^ {xi - x-2)2('^-i) • 

When d = 2> the above computation has already been done in [29] in the context of holographic 
Wilson loops. Moreover, notice that our calculation is essentially the same as that of a two point 
function of a marginal scalar operator (dual to a massless scalar field) in the AdSrf/CFT(i„i setup 
[2]. This suggests that the result (I6.19|) may be interpreted as a two point function of certain 
operators with the conformal dimension d—1 in the dual CFT,^. We naturally expect that such an 
operator is related to the energy stress tensor. It will be an intriguing future problem to work out 
this precisely. 

It is also possible to calculate the entanglement density when we choose the subsystem A to be 

a small perturbation of a round ball B . The corresponding minimal surfaces can be obtained 

by performing the conformal transformation in the AdS space (see |30|): 

,„ x'' + c''(x2 + z2) 
x^ — ^^ 



1 + 2C- X + c2(x2 + Z^) 

z 



z' = ^-^5 ^. (6.20) 

l + 2c-X + c2 x2 + z2 ^ ' 



The result is simply given by 



n(J7i, f)2) = -TT^TT^ = ..„ .-2„ .,„ • ^(STT) > (6-21) 



Oi -O2 



where the unit vectors r2i,2 parameterize the two positions on S and r is the radial coordinate 
orthogonal to S'^~2_ 
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6.4 Entanglement Density and Metric Perturbation in Higher Dimensions 

The main motivation to study the entanglement density in this paper is to understand its dynamical 
properties by using the holographic calculations. Here we would like to study this in the setup of 
AdS4/CFT3. In particular, we will holographically compute the change of entanglement density 
in the presence of small perturbation of the metric compared with the pure AdS result (j6.19p by 
looking at the pure gravity in four dimensions. We again choose ^ to be a small perturbation 
around the half space. 

In the metric background (12. ip with the perturbation (j2.2p at d = 3, the shift of HEE AS a 
for the surface (I6.10p is given by the following expression up to the linear order of /i^j/ and the 
quadratic order of /: 

R"^ f dxdz 



ASa 



AGn J z 



-,2 



hxx , df ^ -hxx + '^hyy (df\'^ ^ h^x + '^hyy (df^"^ 
2 + "^fe+ 4 \dx + 4 \dz 



(6.22) 



The minimal surface / can be found as (]6.16p by setting d = 2>. 

We can obviously divide (f6:22]l into three parts AS a = ASf^ + AS'J^ + AS^^\ where AS^^ is 
defined to be the term which involves i-th power of /. We can find from ASj^ 

5 AS A _ R^ f dxdz ^ ( Az\xi-x) \ 

Saxi)~2GN7Tj z^ ''"''^'[{{x.-xy + z^)^)- ^^•^''^ 

In the above expression, the metric perturbations h^^ are evaluated at any constant values of y 
and t, corresponding to the definition of the subsystem A. 

(2) 

Moreover, we can evaluate the shift of the entanglement density from AS)^ as follows: 

AnA(t,y,xi,X2) 
6^ASa 



5^x1)6^x2) 

K>2 roD 



27r^GN Jo 
where A and B are given by 



dz / dx[A{x,z)hyy{t,y,z) +B{x,z){hxxit,y, z) + hyy{t,y,z))], (6.24:) 



2 (3(x - x i)^ - z^) (3(a; - xa)^ - z^) + 48z^(x - xi){x - X2) 

yi( X, z) — z 



{{x — Xi)^ + z'^) {{x — X2)^ + Z^) 

B{x,z) = z -3^- ,, . ,,3 . (6.25) 



2 (3(x - xxf - z^) (3(x - X2f - z"^) - 16z2(x - xi)(x - xa) 
((x — xi)^ + z^) ((x — xa)^ + z"^^ 



By using the expression (j6.24p we would like to study the relation between the entanglement 
density Aha and energy momentum tensor T^i,. The latter is related to the asymptotic behavior 
of the metric via (|2.8p and (12. 9p . The full metric perturbation is found as (14. 3p by solving the 
Einstein equation. 
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In order to proceed analytically we assume a time-dependent but translationally invariant per- 
turbation of the pure AdS4, which is explicitly expressed as (??). By taking the Fourier transfor- 
mation of An^ with respect to the time t, we finally obtain the following result: 

AnA{uj,y, xi,X2) 

fr\ roc /"CXD 

= -4^ • |a;|"3/2 . / dz dx[Aix,z)Tyy{uj) + B{x,z)(T.^.^iuj)+Tyy{uj))]z^/^J3/2i\^\z), 

V^ Jo J-oo 



e-'l"l/V| • ((^1^1 - 12)Tyy{oj) + 2Tu{oj)) , (6.26) 



where we set / = |xi — X2|. 

We can Fourier transform An^(a;,y,xi,X2) and express the result as a perturbation at time t 
carried by a Green function to time t, obtaining 

AnAit, y, XUX2) = j di (G«(t - i) ■ Tuii) + Gyyit - i) ■ Tyy{t)) (6.27) 



with 



;4 + 6/^T^ - 24r^ _ {f - 4r^ ^ 

(/2+4r2)3 ' ^''^^> ~ 2(/2+4r2)2 



Gyy{T) = -2-—j^.^^^^^;.^^^^, G'it(T) = ,A___(_ (6.28) 



7 Conclusions 



In this paper, we studied the dynamics of holographic entanglement entropy (HEE) induced by a 
metric perturbation around the pure AdS spacetime. Its dynamics is constrained by the perturba- 
tive Einstein equation. First we solved this to express the HEE in terms of the energy stress tensor 
in the dual GET. Next we obtained constraint differential equations (j3.22p and (j3.23p . which are 
satisfied by the perturbation of HEE for the Einstein-Scalar theory in AdSs. We can regard this 
as a counterpart of Einstein equation from the viewpoint of the GFT. 

Moreover, we evaluated the evolution of entanglement entropy in the pure gravity theory on 
AdS4 and AdSs. In particular we found that the increased amount of HEE AS a is given in terms 
of the energy density Tu via a non-local transformation as in the equation (I4.15p . 

It will be an important and interesting future problem to find a universal constraint equation 
for general setups in higher dimensions and to work out how to take into account higher order 
perturbations systemically in these arguments. 

From the GFT viewpoint, our analysis predicts the behavior of entanglement entropy for excited 
states. As we explicitly show for AdSs and AdS4 setups, the entanglement propagates at the speed 
of light and this is dual to the propagation of gravitational waves in the bulk AdS. 
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At the same time our analysis gives a further support of the first law-like relation |25j when we 
choose the subsystem to be a round ball. Our result shows that we get the first law-like relation 
if the spatially inhomogeneous modulation is small enough compared with the subsystem size i.e. 
kl « 1. In this way we expect a certain fundamental mechanism suggested by this robust first 
law-like relation. It will be intriguing to find a field theoretic explanation of this property, especially 
in higher dimensions. 

We also present a higher dimensional extension of entanglement density, which was first intro- 
duced in |16] for two dimensional field theories. We showed that this quantity is non- negative owing 
to the strong subadditivity. We presented its holographic calculation and extract its perturbation 
for weakly excited states. 

Note Added: After the present paper was listed on the arXiv, independent papers [311 [32l [33] 
appeared. In these papers, first law-like relations when the subsystem is a strip have also been 
analyzed. They showed that the variation of entanglement entropy includes other components of 
energy stress tensors other than Tu- 
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